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Equivalence of Bose-Einstein condensation and symmetry breaking 
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Based on a classic paper by Ginibre [Commun. Math. Phys. 8 26 (1968)] it is shown that whenever 
Bogoliubov's approximation, that is, the replacement of ao and ag by complex numbers in the 
Hamiltonian, asymptotically yields the right pressure, it also implies the asymptotic equality of 
|(ao)P/V and {agao) /V in symmetry breaking fields, irrespective of the existence or absence of 
Bose-Einstein condensation. Because the former was proved by Ginibre to hold for absolutely 
integrable superstable pair interactions, the latter is equally valid in this case. Apart from Ginibre's 
work, our proof uses only a simple convexity inequality due to Griffiths. 
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The equivalence of Bose-Einstein condensation (BEC) 
and a spontaneous breakdown of the gauge symmetry 
related to the choice of a global phase factor in creation 
and annihilation operators is an intriguing problem which 
apparently has not yet been solved rigorously. Most of 
the time the equivalence is tacitly assumed; for exam- 
ple, Hohenberg's celebrated result [1] about one- and 
two-dimensional Bose-systems is often considered as the 
proof of the absence of BEC in low dimensions although 
it proves the absence of symmetry breaking. The equiva- 
lence of BEC and symmetry breaking can be summarized 
as the asymptotic equality of \{ao)\'^/V and {aQao)/V in 
the limit of infinite volume, when the thermal or ground 
state averages are taken in the presence of a symmetry 
breaking field. The operator Oq creates a boson in a one- 
particle state distinguished by the fact that its occupa- 
tion may become macroscopic. The problem is that in 
finite volumes the first quantity is always less than the 
second; if equality installs asymptotically, it is the work 
of a large deviation principle resulting in the equivalence 
of certain statistical ensembles. This question, closely 
related to the validity of the Bogoliubov approximation 
(BA) [2] , raised only a limited interest almost exclusively 
among mathematical physicists [3] -[8]. The equivalence 
of the two notions in the above sense was recently shown 
by the present author for some models with a simplified 
interaction [9], without using BA. In that paper it was 
overlooked that the equivalence for a general interacting 
Bose gas was a more or less direct consequence of a pow- 
erful result about BA by Ginibre, obtained almost forty 
years ago [3]. One may say that Ginibre himself over- 
looked this important implication of his work. This note 
is to present the necessary argument. 



making nearly the minimum assumption to guarantee a 
non-pathological thermodynamic behaviour. The pair in- 
teraction (j) is superstable and weakly tempered, i.e. for n 
particles in a box of volume V the total potential energy 



U{ri, . . . ,r„) 



r,; - r. 



)>-bn + an^/V (1) 



for any set of position vectors in the box, and 

cj){v) < c^or-^''+'^ (2) 

for |r| = 7- > i?. Here a > and b are independent of 
the volume, d is the space dimension and e a-nd R 
are positive constants. Ginibre was interested in proving 
the correctness of BA which consists in treating ao and 
flg as complex numbers. Here Oq creates a boson in the 
one-particle state cpo = 1/VV (but tpo could be any other 
one-particle state). Ginibre's interpretation of BA is as 
follows. Let Fq be the single-mode Fock space spanned 



by the product states 



0,1,2,... and let F' be 



the Fock space built on the orthogonal complement of tpQ 
in the one-particle Hilbert space [13]. Let, moreover, Ppi 
be the orthogonal projection onto F' . Then, for an op- 
erator B on F and a complex number C the Bogoliubov 
approximation of B is 



BoiC) = Pf'AcBA*cPf' 



(3) 



where it is understood that Bq (C) acts only on elements 
of F' . The operator 



At 



^Cag-Caa 



g-|C|V2gCaSg-Cao (4) 



applied to a V' in F' creates the coherent state 



In his classic but not widely known paper Ginibre con- 
sidered a Bose gas with a rather general pair interaction 
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tensor-multiplying If tjji , 11)2 are in F', 

{iPi\Bo{C)\ij2) = (Vi ® C). (6) 

In the first part of his paper Ginibre applied the trans- 
formation (3) onto the density matrix W = e~^^ . The 
Hamiltonian of the system is 

H = T + U-nN- pWiao + a*o) (7) 

where T is the kinetic energy, /i is real, A'^ is the parti- 
cle number operator and the amplitude v of the gauge 
symmetry breaking field is chosen to be real. Because 
is unitary, without the projection the transformation 
preserves norm, trace and positivity. Together with the 
projection the trace of a positive operator cannot but 
decrease. So 

Tr Wo(C) = Tr 'AcWA*c < TrW = Z (8) 

where Tr' is the trace in F' and thus 

log Tr Wo(C) <V-HogZ = l3pv{n,iy) ■ (9) 

Ginibre proved that the limit of the pressure pv exists 
and 

lim sup (/3F)"^ log Tr Wo(C) = Km py = pifJ., f) 

V— »00 Q V—>-QO 

(10) 

showing thereby that with the right choice of C, BA for 
the density matrix reproduces the exact pressure in the 
thermodynamic limit. 

Bogoliubov's original idea was to make the replacement 
in the Hamiltonian. Because ao\C) = C\C), writing H 
as a normal-ordered polynomial of creation and annihi- 
lation operators and applying (3), the resulting Ho{C) is 
indeed the same as the outcome of a simple substitution. 
In coherent states particles can enter in contact with each 
other. While (10) holds also for hard-core interactions, 
to make Ho{C) meaningful, hard cores have to be ex- 
cluded. With the additional condition of the absolute 
integrability of </>, Ginibre proved the analog of Eq. (10). 
The inequalities 

Zo{C) = Tr'e-^"°^'^^ < TrWo(C) < Z (11) 

or 

PoiC) = {(3V)-HogZo{C)<pv , (12) 

derived in [3], suggest that, again, po{C) is to be maxi- 
mized. The main result (Theorem 3) of [3] is 

lim sup po{C, fijiy) = lim pv{lJi,v) = p{n,v) (13) 

which is our point of departure. 

The proof of the equivalence of BEG and symmetry 
breaking in the sense discussed in the introduction is as 
follows. 

(i) When a sequence of convex functions /„ converges 



(necessarily to a convex function /), the sequence //^ of 
derivatives also converges apart from possibly a set a zero 
measure. In particular, the inequalities 

f{x-Q) < liminf /^(x-O) < lim sup (cc-^O) < f'{x+0) 

(14) 

due to Griffiths [10] hold true. There is at most a count- 
able infinite set of x such that f'{x — 0) < f'{x + Q); 
otherwise we have equalities in (14). In the first two ap- 
plications below the subscript n will correspond to V. 
(ii) We embed H into a one-parameter family of auxiliary 
Hamiltonians 

H' =T + U-iiN' - fioNo - i^Vv{ao + a*) (15) 

where A^o = «o<^o and N' = N — Nq. It is easily checked 
that the pressure p'y{n, corresponding to H' has 

all the nice properties shown in [3] for py(/i, i^). Es- 
pecially, for fixed p, and v, p'y is a convex increasing 
and real analytic function of /xq (so it is continuous, 
Py{fj,, IJ,,v) = pvil-i,!^)), and for fixed ^ and fio it is a 
convex even and real analytic function of u. These prop- 
erties, apart from the analyticity but including continu- 
ity, are inherited by the (existing) thermodynamic limit 
p'iH, Ho, J^), so p'iii, H, v) = p(/i, v). Also, 

_ f dpy{i^,no,i^) 
VV ^ 2 dv ^ ' 

and 

{No)fj,,^,o,^ ^ dp'y{ii,HQ,v) 

V diio ^ ' 

where the averages are taken with the density matrix 

g-0H I^Ye^^^ . For any fixed /i and for (/io,^^ 7^ 0) 
in a set fi^ of full Lebesgue measure both dp' /dv and 
dp' /duo exist, and by Eq. (14) 

lim {ao)^,^,„./^={l/2)dp'{y.,iJio,y)ldv (18) 

V— >oo 

and 

lim {No)u„u,o,t./V = dp'{n, Ho, i')/dno ■ (19) 

V ^oo 

With the possible exception of a countable number of 
values of u Eq. (18) holds also for hq = fj, . This may 
not be true for Eq. (19): there is an abstract possibility 
that for a positive-measure set of v the /ig-derivative of 
p' does not exist at /xo = M- This, however, would have 
no practical importance because for any choice of (/i, v), 
dp' /diiQ exists for almost every fiQ and. hence, for fiQ 
arbitrarily close to fi. We return to this point later, 
(in) Apply BA to H' to obtain i?^(C), Z^(C) and p(,(C). 
Theorem 3 of [3] extends without any further ado to this 
case, resulting p'o{C) < p'y and 

lim sup p'o{C) =p' (20) 

V—t-OO (J 
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for all /i, /ig and v. With the choice v ^ Q real the 
maximum is attained for C real, Cv > 0. Within this set 
C, p'q reads 

p'o{C,n,fio,iy)^p'o{C,fi,0,0)+fio^+2\i.\^ . (21) 

For fixed C and /x this is a convex function of both /xq 
and v; 

SUpPo(C) = Po(Cinax(M, Mo, i^), M, Mo, z^) (22) 

c 

is the upper envelope of {po(C, /i, /ig, J^)|C G C} and is, 
therefore, convex in /xq and The finite- volume pressure 
Po(C, /X, 0, 0) = Pq{—C, /X, 0, 0) is a real analytic function 
of C. So p^(C, 1^, 0, 0) « p^(0, /X, 0, 0) + aC^ close to C = 
0. Thus, if 1/ 7^ 0, Cmax 7^ either, and dp'o/dC = at 
C = Cmax- This implies 

{lJ,,Un,l/),U,Llr),l/) 

C'max 

i = 2 sgn U ■ 



and 



^Po(<^max(M>/^0,J'),M>MO,J') _ C*, 



2 

max 



(23) 
(24) 



a/xo 1^ 

Now we use (14) a second time, with /„ = 
Po(Cmax,M.Mo,J^) and / = p\iJL,iJiQ,v). For any ji and 
any (/xo,2^) e 0^:^ 

/xo,J/)= lim |Ci„ax(M,Mo,J^)|/\/^ (25) 

V — »c» 

exists, and 

hm = sgni/c^3^ = sgnv lim 

V — »oo \/ Y V — >oo 



(-/Vo)p,po,t/ 



(26) 

(iv) The above result can bo improved. When 
dp{^^v) / dp exists, we can set /xq = /i in Eqs. (18), (23) 
and (25) and obtain the first half of Eq. (26), 

ldp{n,v) (ao)^,L/ , ^ 

2^^ = ^^Too^=^Sn^w(M,^). (27) 

Here Cmax(/^-'^) = c[,jjj^(^, ^, i/) which is clear from 
Po{C,IJ,,v) = p'q{C, ^, fi,iy). Now fixing any /i, for all 
p outside a possible (/x-dependent) set of zero measure 
there exist both the derivative dp{iJ,, v) /dv and a se- 
quence /in i such that (/i„,j/) G f2^. For such 
a v and /i„, p'{^, ^n,i^) ^ p{l^,v) by continuity and 
dp'{fi, fin,i^)/di' dp{ii,v)/ dv by Griffiths' Lemma 
(14). Combining Eqs. (18), (26) and (27) one obtains 
that for almost all {^i, p ^ 0) 



lim 



- sgn i/ Cmax (A«,J^) 



= sgnu lim lim .liMjh!^, (28) 

n— *oo y— »oo V K 

Equation (28) is nearly the final result except that in 
the right member the two limits should be interchange- 
able. This, however, is easily seen. Because of the 
Schwarz inequality. 



< lim 



V 



(29) 



Suppose that there is a strict inequality here. Note that 
(-^o)/i,/io,i^ is an increasing function of /Xq. Then for some 
S > and for all F > V5 and all /(xq > /U 



l^^o)'*.''! +^ < ^/IM/^ < ^/IM^ii^foi^ . (30) 



V 



V 



The contradiction is obvious because /U„ J, /x in Eq. (28). 
Our final result is, thus, that for almost all (/x, v ^0) 



lim 



sgn IV Cmax (M> 2^) 



: sen u lim 



{No) 



V 

(31) 

Equation (31) proves that an optimized c-number substi- 
tution provides the right answer at least for the asymp- 
totic value of {aQ)/VV and {aQaa)/V. It docs not tell 
anything about these values (although for i^ ^ they 
are surely of order 1, as in [9]), especially, about their 
limit when v tends to zero. So the equation does not 
prove that there is a spontaneous symmetry breaking, 
but it does prove that symmetry breaking occurs simul- 
taneously with Bosc-Einstein condensation. 

In approximate calculations one often starts with the 
canonical transformation oq = 60 + (do) where, by def- 
inition, the operator has zero mean; sec e.g. [11]. 
The approximations made afterwards arc based on the 
hypothesis that the fluctuations of 60 arc negligible or 
small compared with (a'^aa). Equation (31) justifies this 
hypothesis by proving (6S^o)m,i^/^ — > as F — > 00. 

The functional form of p' is 

p'{fl, HQ, V) = .9(Cmax(M, Mo, v), M) 

+/«ociax(M>Mo,J^) + 2|i/| {h,Hq,v) . (32) 
Alternately, p' has to satisfy the equation 



dp'{ii, Ho, v) If dp'{n, Mo, v 



9/xo 



dv 



(33) 



obtained by reading the two ends of Eq. (26). Apart 
from the fact that dp' / dv tends to dpjdv, as /xo tends to 
/tx, the Bose gas described by the extended Hamiltonian 
H' may have its own interest. In any case, if /xq is close 
to /i, the behaviour of this system is not expected to 
significantly differ from that of the original one. BEG or a 
spontaneous symmetry breaking is described by solutions 
of (33) that arc convex in all the three variables and have 
a discontinuous i/-derivative at = 0. 

A different solution to the problem treated in this Let- 
ter can be found in Ref. [12]. 

This work was partially supported by OTKA Grants 
T 042914 and T 046129. 
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